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Abstract 

We have experimentally checked the Jarzynski equality and the Crooks relation U[2||2] on 
the thermal fluctuations of a macroscopic mechanical oscillator in contact with a heat reservoir. 
We found that, independently of the time scale and amplitude of the driving force, both 
relations are satisfled. These results give credit, at least in the case of Gaussian fluctuations, 
to the use of these relations in biological and chemical systems to estimate the free energy 
difference AF between two equilibrium states 21 El- An alternative method to estimate AF 
in isothermal process is proposed too. 

1 Introduction 

Many systems in Nature and in technological applications work out of equilibrium. However, a 
precise estimation of the free energy difference AF between two equilibrium states A and B of these 
systems is extremely useful to increase our knowledge of the underlying physical processes which 
control their dynamical behaviour. It is well known that AF can be estimated by perturbing the 
system with an external parameter A and by measuring the work W done to drive the system from 
A to B. However this method gives in general an overestimation of AF because W > AF, where 
the equality holds if and only if the perturbation is infinitely slow. In other words the path 7 to go 
from AtoB has to be a reversible one. In many systems, because of unavoidable experimental and 
environmental constraints, the path 7 is not a reversible one, that is the system cannot be driven 
from A to B in a time much longer than its relaxation time. This may happen for example in all 
the systems where thermal fluctuations cannot be neglected and the external power injected into 
the system is comparable to the thermal energy. In 1997 fl] Jarzynski derived an equality which 
relates the free energy difference of a system in contact with a heat reservoir to the probability 
distribution function (pdf ) of the work performed on the system to drive it from A to B along any 
path 7 in the system parameter space. Speciflcally, when A is varied from time t = to t — tg, 
Jarzynski deflnes for one realization of the "switching process" from A to S the work performed 
on the system as 

dH4z{t)] 



I X—^dt, (1) 
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where z denotes the phase-space point of the system and li\ its A-parametrized Hamiltonian (see 
also j^). One can consider an ensemble of realizations of this "switching process" with initial 
conditions all starting in the same initial equilibrium state. Then W may be computed for each 
trajectory in the ensemble. The Jarzynski equality (JR) states that 

AF=-^log(exp[-/3M/]), (2) 
P 

where (•) denotes the ensemble average, /9~^ = ksT with fc^ the Boltzmann constant and T the 
temperature. In other words (exp [— /3Wdiss]) — 1, since we can always write W — AF+Wdiss where 
Wdiss is the dissipated work. Thus it is easy to see that there must exist some paths 7 such that 
W^diss < 0. Moreover, the inequality (expx) > exp (a:) allows us to recover the second principle, 
namely (Wdiss) > 0, i.e. (W) > AF. From an experimental point of view the JE is quite useful 
because there is no restriction on the choice of the path 7 and it overcomes the above mentioned 
experimental difficulties. Numerous derivations of the JE has been produced [3 |H1 El EH E] > but it 
seems from the recent criticisms of Cohen and Mauzerall ^2] that this result is still under debate. 
Without willing to enter in this letter into the theoretical debate, we think that it is important 
to experimentally check the JE on a very simple and controlled system in order to safely use it 
in more complex cases as the biological and chemical ones, where it is much more difficult to 
verify the results with other methods. For this reason we experimentally probe a model system: a 
macroscopic mechanical oscillator driven out of equilibrium, between two equilibrium states A and 
B, by a small external force. We show that the JE is experimentally accessible and valid, and does 
not depend on the oscillator's damping, on the driving force's switching rate and on its amplitude. 
In our experiment we can also check the Crooks relation (CR) which is somehow related to the JE 
and which gives useful and complementary information on the dissipated work. Crooks considers 
the forward work Wf to drive the system from ^4 to i? and the backward work Wb to drive it from 
B to A. If the work pdfs during the forward and backward processes are Pt(W) and Pb(VF), one 
has [lEI 

^^^^ = exp {(i\W - AF]) - exp \m^..\ (3) 

A simple calculation from Eq. lO leads to Eq. iQ. However, from an experimental point of view 
this relation is extremely useful because one immediately sees that the crossing point of the two 
pdfs, that is the point where Pf(VF) = Pb(— W), is precisely AF. Thus one has another mean 
to check the computed free energy by looking at the pdfs crossing point AFx . Let us examine 

in some detail the Gaussian case: P(W^) oc exp ^— ^ leads to AF — (W) — i.e. 

(Wdiss) = ^^Y' ^ ^- Furthermore, it is easy to see from Eq. ^ that if Pf(W^) and Pb(— VF) are 
Gaussian, then AF = and /3cr^ = (VF)f + (VF)b = 2 (Wdiss). Thus in the case of 

Gaussian statistics AF and Wdiss can be computed by using just the mean values and the variance 
of the work W. 

Before describing the experiment, we want to discuss several important points. The first is the 
definition of the work given in Eq. iQ, which is not the classical one. Let us consider, for example, 
that A is a mechanical torque M appHed to a mechanical system S, and —dH\/d\ the associated 
angular displacement 9. Then, from Eq. (QJ, one has 



W = - MOdt 







M9 



where W"^ ^ - I MOdt (4) 
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Figure 1: (i) Schematic drawing of the oscillator; (ii) FDT check in oil: the circles (o) correspond 

to the direct measurement of the noise and the dashed curve is ^ (|^| ) computed by inserting the 
measured response function in the FDT, Eq. 10 



is the classical work. Thus W and W^^ are related but they are not exactly the same and we 
will show that this makes an important difference in the fluctuations of these two quantities. The 
second point concerns the Ai^ computed by the JE in the case of a driven system, composed by S 
plus the external driving. The total free energy difference is Ai^ — AFq — $ where AFq is the free 



energy of S and $ — 



MB 



the energy difference of the forcing. The JE computes the Ai^ of the 

A 

driven system and not that of the system alone which is AFq . This is an important observation in 
view of all applications where an external parameter is added to S in order to measure AFq [S]. 
Finally we point out that, in an isothermal process, A_Fo can be easily computed, without using 
the JE and the CR, if W^^ is Gaussian distributed with variance tT^c . Indeed the crossing point 
of the two Gaussian pdfs Pf(H/'=') and Pb(-1^'=') is 



which by deflnition is just -AfJ), i-e. = -AFq. Furthermore {W''^){ + = -2 (PKiis 



by deflnition, but in this case 2 (VFdiss) 7^ P<^1 



2 Experimental setup 

To study the JE and the CR we measure the out-of-equiHbrium fluctuations of a macroscopic 
mechanical torsion pendulum made of a brass wire, whose damping is given either by the vis- 
coelasticity of the torsion wire or by the viscosity of a surrounding fluid. This system is enclosed 
in a cell which can be fllled with a viscous fluid, which acts as a heat bath. A brass wire of 
length 10 mm, width 0.75 mm, thickness 50 /im, mass 5.91 x lO^'^ g, is clamped at both ends, 
hence its elastic torsional stiffness is C = 7.50 x 10^'* Nmrad~^. A small mirror of effective mass 
4.02 X 10^^ g, length 2.25 mm, width 6i = 7 mm, thickness ai — 1.04 mm, is glued in the middle 
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Figure 2: Oscillator immersed in oil [case (a)]: (i) Applied external torque, (ii) Induced angular 
displacement, (iii) its psd, (iv) its pdf, (v) Injected power computed from the Jarzynski definition 
W = —MO, (vi) Injected power computed from the standard definition VV^^ = —MO 



of the wire, see Fig.nji), so that the moment of inertia of the wire plus the mirror in vacuum 
is / = 1.79 X 10~^° kgm^ (whose main contribution comes from the mirror). Thus the resonant 
frequency of the pendulum in vacuum is /o — 326.25 Hz. When the cell is filled with a viscous 
fiuid, the total moment of inertia is leS = / + /fluid, where /fluid is the extra moment of inertia given 
by the fiuid displaced by the mirror p^l- Specifically, for the oil used in the experiment (which 
is a mineral oil of optical indice n — 1.65, viscosity u — 121.3 mPas and density p — 0.9pwatcr 
at T = 21.3 °C) the resonant frequency becomes /o = 213 Hz. To apply an external torque M 
to the torsion pendulum, a small electric coil connected to the brass wire is glued in the back of 
the mirror. Two fixed magnets on the cell facing each other with opposite poles generate a static 
magnetic field. We apply a torque by varying a very small current J fiowing through the electric 
coil, hence M cx J. The measurement of the angular displacement of the mirror is done using a 
Nomarski interferometer JHEj whose noise is about 6.25 x 10^^^ rad/\/Hz, which is two orders 
of magnitude smaller than the oscillator thermal fluctuations. A window on the cell allows the 
laser beams to go inside and outside. Much care has been taken in order to isolate the apparatus 
from the external mechanical and acoustic noise, see QS] for details. 

The motion of the torsion pendulum can be assimilated to that of a harmonic oscillator damped 
by the viscoelasticity of the torsion wire and the viscosity of the surrounding fluid, whose motion 
equation reads in the temporal domain 

hsO+ I G{t-t')0{t')dt' + C0 ^ M, (6) 

where G is the memory kernel. In Fourier space (in the frequency range of our interest) this 
equation takes the simple form [—IcS uj'^ + C]0 = M, where • denotes the Fourier transform and 
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Figure 3: Case a): (i) Pf{W) and Pb(-VF), (ii) PfiW^) and Pfi-W^); Case c): (iii) Pi{W) and 
Ph{—W), (iv) Pf(M^'^') and P{{—W^^) (experimental forward and backward pdfs are represented 
by o and □ respectively, whereas the continuous lines are Gaussian fits) 



C — C + i[C'{ + wCj ] is the complex frequency-dependent elastic stiffness of the system. C'{ and 
C2 are the viscoelastic and viscous components of the damping term. The response function of 
the system x — 0/ M can be measured by applying a torque with a white spectrum. When M = 0, 
the amplitude of the thermal vibrations of the oscillator is related to its response function via 
the fiuctuation-dissipation theorem (FDT) 0. Therefore, the thermal fiuctuation power spectral 
density (psd) of the torsion pendulum reads for positive frequencies 

(1^1^) = l^lmx = 91±^ . (7) 

We plot in Fig.njii) the measured thermal square root psd of the oscillator. The measured noise 
spectrum [circles in Fig.QIii)] is compared with the one estimated [dotted line in Fig.QIii)] by 
inserting the measured x in the FDT, Eq. l(£||. The two measurements are in perfect agreement 
and obviously the FDT is fully satisfied because the system is at equihbrium in the state A where 
M = (see below). Although this result is expected, this test is very useful to show that the 
experimental apparatus can measure with a good accuracy and resolution the thermal noise of the 
macroscopic pendulum. 
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3 Experimental results 



Now we drive the oscillator out of equilibrium between two states A (where M — 0) and B 
(where M = Mmax = const ^ 0). The path 7 may be changed by modifying the time evolution 
of M between A and B. We have chosen either linear ramps with different rising times t, see 
Fig.|2li), or half-sinusoids with half-period t. In the specific case of our harmonic oscillator, 
as the temperature is the same in states A and B, the free energy difference of the oscillator 



alone is AFq = AU ^ 



whereas AF — AFq 



i.e. for an harmonic 



2C 

potential AF = —AFq. Let us first consider the situation where the cell is filled with oil. The 
oscillator's relaxation time is t^-cUx — 23.5 ms. We apply a torque which is a sequence of linear 
increasing /decreasing ramps and plateaux, as represented in Fig.EIi). We chose different values of 
the amplitude of the torque M [11.9, 6.1, 4.2 and 1.2 pNm] and of the rising time r [199.5, 20.2, 
65.6, 99.6 ms, respectively], as indicated in Tabled [cases a). . .e)]. Thus we can probe either the 
reversible (or quasi-static) paths (r » Treiax) or the irreversible ones (r < Tioiax)- We tune the 
duration of the plateaux (which is at least 4 Treiax) so that the system always reaches equilibrium 
in the middle of each of them, which defines the equilibrium states A and B. We see in Fig.|2Iii), 
where the angular displacement 9 is plotted as a function of time [case a)], that the response of the 
oscillator to the applied torque is comparable to the thermal noise spectrum. The psd of 9 is shown 
in Fig.[2liii). Comparing this measure with the FDT prediction obtained in Fig.QJii), one observes 
that the driver does not affect the thermal noise spectrum which remains equal to the equilibrium 
one. Moreover we plot in Fig.[2liv) the pdf of the driven displacement 9 shown on Fig.|2lii), which 
is, roughly speaking, the superposition of two Gaussian pdfs. From the measure of M and 9, the 
power injected into the system W can be computed from the definition given in Eq. 1^, that in 
this case is = —M9. Its time evolution, shown in Fig.|2lv), is quite different from that of the 
classical power W'^^ — —M9, whose time evolution is plotted on Fig.[2tvi): W is non-zero only 
for M and vice- versa W"^^ ^ only for M ^ 0. From one time series of W we can compute 
from Eq. (QJ the forward and the backward works, W[ and Wb, corresponding to the paths A ^ B 
and B ^ A, respectively. We also do the same for the classical work. We then compute their 
respective pdfs Pf(M^) and Pb(— VF). These are plotted on Figs.|2li,iv) where the bullets are the 
experimental data and the continuous fines their fitted Gaussian pdfs. In Figs.El the pdfs of W and 
W^^ cross in the case a) at PW ~ —23.5, and in the case c) PW ~ —6.1. These values correspond 

to the AF = — = - AFq- We find that this result is true independently of the ratio r/riciax 
and of the maximum amplitude of |M|, Mmax- This can be seen on Tabled where the computed 

AU = is in good agreement with the values obtained by the crossing points of the forward 

and backward pdfs, that is AFx for P(W) and —AW^^ for P(VF'^'). Finally inserting the values 
of Wf and Wh in Eq. ^ we directly compute AF[ and AFb from the JE. As can be seen in Table 
n] the values of AFq obtained from the JE, that is either — (AFf -I- <I>) or — (AFb -|- <&), agree within 
experimental errors with the computed AU. Indeed JE works well either in the foreseeable case 
a) where t ^ Trdax or in the critical case b) where r < Treiax- The other case we have studied 
is a very pathological one. Specificafiy, the oscifiator is in vacuum and has a resonant frequency 
fo — 353 Hz and a relaxation time Trdax = 666.7 ms. We applied a sinusoidal torque whose 
amplitude is either 5.9 x 10~^^ or 9.4 x 10"^^ Nm [cases f) and g) in Tabled respectively]. Half 
a period of the sinusoid is t = 49.5 ms, much smaller than the relaxation time, so that we never 
let the system equilibrate. However, we define the states A and B as the maxima and minima 
of the driver. Surprisingly, despite of the pathological definition of the equifibrium states A and 
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Table 1: Free energies of cases a)...g) defined in the text (the values of Mmax are in pN m). 
AJ7 = is the computed expected value 

-B, the pdfs are Gaussian and the JE is satisfied as indicated in Table Moreover, this happens 
independently of Mmax and of the critical value of the ratio r/rrciax ^ 1- Finally, we indicated in 
Table the value AFq which is the free energy computed from the JE if one considers the "loop 
process" from A to A (the same can be done from B to B and the results are quantitatively the 
same). In principle this value should be zero, but in fact it is not since we have about 5% error in 
the caHbration of the torque M. 

4 Conclusion 

Our results clearly demonstrate the validity and the robustness of the JE in an isothermal process, 
at least when the work fluctuations are Gaussian and when the harmonic approximation is relevant 
for the system. The more accurate and reliable AF estimator is given by the crossing points AFx 
and W^, because they are less sensitive to extreme fluctuations which may perturb the convergence 
of the JE. We have also shown that, in the case of Gaussian fluctuations, remains an excellent 
estimator even in cases where the JE and the CR could not hold. Unfortunately our results do not 
fully throw light on the debate between Cohen, Mauzerall and Jarzynski j2lE| since the work pdfs 
are Gaussian. Recently, Ritort and coworkers have used the JE to estimate AF in an experiment 
of RNA stretching where the oscillator's coupling is non-linear and the work fluctuations are non- 
Gaussian It would be interesting to check these results on a more simple and controlled 
system. We are currently working on the experimental realization of such a non-linear coupling, 
for which AF ^ —AFq. Finally we want to stress that our results, although limited to the Gaussian 
case, show that it is possible to measure tiny fluctuations of work in a macroscopic systems. As 
consequence they open a lot of perspective to use the JE, the CR and the recent theorems on 
dissipated work (see for example to characterize the slow relaxation towards equilibrium in 
more complex systems, for example aging materials such as glasses or gels |2(1] . 
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